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'  ABSTRACT 

Data  collected  during  a  continuous,  surface  release,  point  source 

/ 

tracer  experiment  off  the  California  coast;,  is  analyzed.  The  effects  of 


high  speed  data  collection  frtrn  an  airborne  platform  are  removed  by  inverse 

transformation  using  fefie  collecting  instrument's  transfer  function, 
in  frequency  space.  jiShe  tracer  plane  is  characterized  by  a  variety  of 

parameters,  including  the  ccxxventicnal  hourly  averaged  sigma-y  and  sigma-z 
values  widely  used  in  Gaussian  plume  dispersion  farnulae.  Gaseous  dispersion 
is  parameterized  for  the  over-water  case  by  classifying  the  tracer  results 
by  stability  in  a  Pasquill-Gif ford  equivalent  scheme,  and  analytically 
describing  horizontal  and  vertical  plume  growth  as  a  function  of  plume 
travel  distance.  Several  other  over-water  data  sets  are  used  in  this 
parameter! zat ion  -  Ccnpariscns  are  made  to  the  over  land  case. 
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INTRODUCTION 


^^The  Minerals  Management  Serviced  formerly  the  Bureau  of  Land 
Management) ^sponsored  a  3eries  of  four  atmospheric  tracer  experi¬ 
ments  at  California  coastal  locations  over  a  two-year  span,  1980- 
1982.  These  experiments  were  designed  to  assess  air  pollution 
impact  from  proposed  oil  exploration  and  drilling  activities  along 
the  continental  shelf.  Two  experiments  (winter  and  summer)  at 
each  of  two  sites  (open  coast  and  Santa  Barbara  Channel)  were 
funded  in  order  to  investigate  air  quality  impact  under  a  range  of 
meteorological  conditions  and  sites. ^  A  brief  summary  of  these 
experiments  and  references  is  supplied  in  Table  1. 

-—zy  The  basic  designs  of  all  four  experiments  were  similar.  A 
tracer  gas,  100%  SF6,  was  continuously  released  from  a  stationary, 
sea  surface  platform  located,  for  the  majority  of  the  experiments, 
approximately  3  miles  from  shore.  During  parts  of  the  last 


expenmi 


iment ,  the 

-Mo  f  UJ 


platform  was  moved  to  distances  up  to  5  miles  from 


shore.  A  variety  of  meteor logical  parameters  were  continuously 
monitored  at  various  locations.  Tracer  gas  concentrations  were 


measured  by  a  variety  of  methods  at  positions  downwind  of  the 
release  platform,  with  the  majority  made  on  or  near  the  shore 
since  the  purpose  was  to  assess  potential  on-shore  air  pollution 


impact.  Experiments  were  limited  to  daytime  periods  of  on-shore 
flow.  Meteorological  measurements,  however,  were  not  restricted 
to  those  time  periods.  This  report  utilizes  a  subset  of  the  data 
base  collected  during  the  fourth  experiments:  offshore,  aircraft, 


Table  1 


Central  California  Coastal  Air 
Quality  Studies,  1981-1982 
(Sponsored  by  Mineral  Management  Service) 


DATE  LOCATION  FINAL  REPORTS  AVAILABLE 


Sep  80  Santa  Barbara 
Channel  Area 

Jan  81  M 

Dec  81  Pismo  Beach  Area 


Aerovironment,  Inc. 


Stanford  Research 
Institute 


Jun  82 


Stanford  Research 
Institute  and  Naval 
Postgraduate  School 


♦Other  reports  available. 


REF.* 


Zanett  i 
et  al.  81 


Dabberdt 
et  al.  83 

Dabberdt 
et  al.  83 
Other 


+NPS  work  was  sponsored  by  both  the  Minerals  Management  Service 
and  the  NPS  Foundation  Research  Program. 


continuous  gas  analyzer  measurements . 

The  intention  of  this  report  is  to  characterize  over-water 
diffusion  from  a  continuous,  near-surface,  point-source  release 
based  upon  these  measurements.  This  report  is  built  upon  the 
meteorological  results  of  Schacher  et  al.  (1982)  and  a 
preliminary  tracer  gas  and  ranging  results  of  Schacher  et  al . 
(1983).  For  a  detailed  description  of  the  measurement  techniques 
and  data  description,  the  reader  is  referenced  toward  these 
reports. 

OUTLINE 

A  report  flow  chart  is  provided  in  Figure  1 .  This  document 
is  organized  into  two  chapters  with  distinctly  different  designs. 
Chapter  1  contains  technical  procedures  and  data  used  in  the 
piece-wise  analysis  of  the  data  set.  The  second  chapter  presents 
one-hour  average  plume  dispersion  parameters,  ay  and  az,  as  a 
function  of  the  well-known  Pasquill-Gif ford  stability  categories 
adapted  for  overwater  use.  Some  additional  data  from  other 
experiments  supplement  our  data  set  to  produce  a  more  general 
parameterization . 

Readers  interested  primarily  in  plume  dispersion  over  water 
are  advised  to  skip  most  of  Chapter  1,  concentrating  mainly  on 
Step  6  and  Chapter  2.  Those  readers  interested  in  the  particular 
techniques  used  in  the  analysis  of  tracer  data  obtained  from  a 
high-speed  platform  may  be  more  interested  in  Step  2.  In 
addition  to  one-hour  standard  <jy  values,  a  wide  variety  of 


Figure  1. 


analytical  exp¬ 
ressions  for 
Cy  and  Ox  as 
a  function  of 
stability  for 
surface  release 
and  ranges  to 


CHAPTER  1  -  DATA  ANALYSIS 
Step  1  -  Organization 

The  following  data  analysis  was  performed  in  a  step-wise 
fashion,  with  the  complete  data  set  stored  and  cataloged  in  the 
Naval  Postgraduate  School's  IBM  3033  mass  storage  system  and 
9-track  tape  at  the  completion  of  each  step.  Performing  the 
overall  analysis  in  six  separate  steps  allowed  for  manual 
interrogation  of  the  data  set  at  each  fundamental  level  and  will 
allow  for  easy  and  flexible  re-analysis  of  the  data  set  in  the 
future. 

The  analysis  starts  with  SF6  concentrations , represented  as 
digital  voltage  output  from  a  continuous  gas  analyzer  for  single 
passes  through  the  plume.  Aircraft  position  was  recorded  from 
dual  miniranger  signals,  resulting  in  paired  position/  concentra¬ 
tion  data.  Each  plume  transect  wis  chosen  to  start  where  the 
analyzer  first  sensed  SF6  along  the  flight  path. 

The  original  data  set  consists  of  seven  files  with  one 
experimental  day  per  file.  Each  file  contains  a  different  number 
of  passes.  Each  pass  starts  a  new  (2048-byte)  data  block?  the 
number  of  blocks  needed  depends  on  the  pass  length.  Records  are 
of  variable  length. 

This  data  set  was  written  into  mass  storage  on  the  IBM  3033. 
The  type  of  mass  storage  file  used  for  this  analysis  is  called  a 
"partitioned  data  set".  This  data  set  consists  of  a  number  of 


user-specified  "members". 


Each  member  can  be  accessed  interac- 


tively  or  via  program  control.  If  the  members  are  less  than  5000 
lines,  they  can  be  edited  interactively  by  the  user.  This  was 
desirable;  therefore,  care  was  taken  to  keep  each  member  under 
this  limit.  Also,  members  must  consist  of  80-character  records. 
Therefore,  the  initial  records  became  unsynchronized  with  the  mass 
storage  records  after  the  transfer. 

At  this  point,  a  simple  program  named  FORMAT  converted  the 
variable  length  record  format  to  a  fixed  length  format.  The 
output  was  interrogated  and  calibration  passes+  removed. 
Calibration  factors  derived  were  added  to  the  header  of  each  pass . 
These  data  were  written  to  a  partitioned  data  set  named  AIR2, 
residing  on  the  Environmental  Physics  Group's  private  mass 
partitioned  storage  volume.  All  data  set  member  names,  format, 
and  content  will  be  presented  in  tabular  form  later  in  this 
report. 

Next,  the  data  set  AIR2  was  transformed  to  AIR3  by  the 
program  REDUCE.  This  program  performed  3  vital  functions. 

First,  it  converted  raw  voltages  (corrected  for  background  SF6 
concentration)  to  parts  per  trillion  (PPT)  concentration  via  the 
calibration  factors  mentioned  above  and  experimentally-derived 
calibration  formulae.  The  calibration  factors  were  periodically 
measured  during  the  experiment.  The  conversion  formulae  account 
for  instrument  non-linearity  at  high  concentrations .  The 
equations  are: 

41  During  a  calibration  pass,  the  instrument  was  purged  with  a 
"span"  gas  of  known  concentration  in  the  instrument's  linear 
region  to  obtain  calibration  factors. 


c 

where  V  is  voltage  normalized  to  laboratory  conditions; 
VQ  is  output  voltage  from  the  analyzer; 

V  is  baseline  (background)  voltage; 

D 

C  is  the  calibration  factor  determined  during  the 
experiment  (See  Table  2). 


SF6  =  5340V  [V  <  1.345]  (2a) 
SF6  *  exp  (1.160V2  -  2.455  V  +  10.122)  Cl. 345  <  V  <  1.687]  (2b) 
SF6  =  exp  (1.461V  +  6.823)  [1.687  <  V  <  2.053]  (2c) 
SF6  =  exp  (4.252  V2  -  16.780  V  +  26.369)  [V  >  2.053]  (2d) 

SF6  is  concentration  in  parts  per  trillion. 
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Calibration  Factors 


6/21/82 

BEGIN 

1640 

6/22/82 

BEGIN 

1720 

6/24/82 

BEGIN 

6/25/82 

Begin 

1300 

1345 

1440 

1520 

6/27/82 

BEGIN 

1720 

6/28/82 

BEGIN 

6/29/82 

BEGIN 

The  second  vital  function  performed  by  “REDUCE"  was  to 
determine  plume  transect  Cartesian  coordinates.  This  was 
accomplished,  in  most  cases,  with  the  mini-ranger  data.  Three 
scenarios  existed,  depending  on  mini-ranger  performance  for  a 
given  pass.  When  both  mini-ranged  distances  were  available, 
polyncminal  fits  were  performed  to  eliminate  data  “jitter"  and 
simple  triangulation  used  to  determine  plume  coordinates.  When 
one,  or  both,  mini-ranger  signals  were  intermittent,  regression 
analysis  was  used  where  possible,  to  fill  in  the  "gaps".  When  one 
or  both  mini-ranger  signals  were  missing,  coordinate  determination 
was  postponed  for  later  analysis.  An  in-depth  discussion  of  the 
above  process  design  is  given  in  Schacher  et  al.  (1983). 

The  sampling  grid  coordinate  system  is  shown  in  Figure  2. 

The  mini-ranger  transmitters  were  located  on  the  beach,  and  north 
and  south  buoys  were  located  so  as  to  aid  aircraft  nagivation. 
Their  grid  map  locations,  along  with  the  variable  ship  locations 
are  given  in  Table  3. 


Figure  2. 


NPS  sf6  tracer  study 

SAMPLING  GRID 
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Table  3 


Grid  Map  Locations 


Reference  Grid  Coordinates 

(See  Figure  2)  _ (meters) 

X  Y 


N .  Bouy 
S .  Bouy 
S'-iip  6/21/82 

”  6/22/82 

"  6/24/82 

"  6/25/82 

"  6/27/82 

M  6/28/82 

'*  6/29/82 

Estuary 

N.  Ref.  Station 
S.  Ref.  Station 


5926 

11140 

8114 

7550 

4055 

10200 

4945 

6369 

4103 

8628 

4111 

8601 

399 

11090 

581 

12493 

1120 

10459 

8896 

15430 

10000 

10000 

10000 

7050 
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INDIVIDUAL  LINE  KEY:  column  one  is  the  code 

P:  date,  time,  pass  number,  data  quality  index,  altitude 
B:  total  elapsed  time  (sec),  total  number  of  points, 

mini-ranger  0  number  of  points,  mini-ranger  1  number 
of  points,  number  of  bad  points,  calibration  factor 
D:  mini-ranger  code,  distance  from  mini-ranger  (m) , 
analvzer  output  (volts) 


Figure  3b. 

MASS  STORAGE  DATA  SET  "Am"; 
ABBREVIATED  SAMPLE  OF  ANALYSIS  STEP  1  (continued) 
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header  information  (self  explanatory) 
mini-ranger  statistics 

R/V  ACANIA  position  in  grid  (see  Figure  2) 
data  -time  in  seconds 

-width  is  distance  from  first  non-zero 
concentration  in  meters 
-concentration  in  PPT 
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Step  2  -  Data  Transformation 

This  step  accounts  for  instrument  effects  on  the  data.  If  the 
data  are  perceived  as  a  time  series,  and  we  treat  the  instrument  as 
a  first-order  linear  system*,  then, 

«v 

flX  'V'V 

. — -  =  AX  +  BU  (3) 

dt 

where  X  is  a  one-dimensional  matrix  of  state  variables; 

0  is  a  matrix  inputs; 

A,  B  are  square  matrices  of  coefficients;  and 
t  is  the  independent  variable. 

In  general,  the  system  output  can  be  represented  as  a  linear 
combination  of  the  state  variables  and  the  inputs.  In  this  speci¬ 
fic  case,  the  input  is  the  true  SF6  concentration,  and  the  output 
of  interest  Jls  a  state  variable;  the  measured  SF6  concentration. 
Also,  this  case  is  concerned  with  only  one  state  variable;  there¬ 
fore,  matrix  expressions  will  be  dropped.  A  solution  can  be 
expressed  as  the  convolution  of  the  input  waveform  and  a  function 
called  the  unit  impulse  response  of  the  system. 


x(t)  =  _/+*h(t  -  t)u(t)dt  (4) 

vdiere  x(t),  u(t)  are  singular  state  and  input  variables, 

respectively; 

h(t)  is  the  unit  impulse  response. 

*In  a  second-order  system,  a  second  state  variable  would  simply  be 
the  derivative  of  the  first  variable. 


Convolutions  are  rather  difficult  to  perform  on  digital 
machines;  therefore,  we  use  the  convolution  theorem,  which  states 
that  convolution  in  the  time  domain  is  analagous  to  multiplication 
in  the  frequency  domain. 

x(t)  -  F_1[X(f)]  =  F-1[H(f)U{f)]  (5) 

where  X(f)  is  the  Fourier  transform  of  x(t' 

U(f)  is  the  Fourier  transform  of  ufr'' 

H(f)  is  the  Fourier  transform  of  h(t),  or  the  "transfer 
function" 

F-1  refers  to  the  inverse  Fourier  transform. 

Finally,  since  the  system  input  is  the  desired  quantity. 
Equation  (5)  is  inverted,  yielding: 

u ( t )  -  F-1CX(f)/H(f)]  (6) 

A  graphical  example  can  serve  as  a  "proof"  of  this  concept, 
often  referred  to  as  the  transfer  function  approach  (see  Figure 
4).  Let  the  unit  impulse  response  be  the  unit  impulse.  The 
impulse  transforms  to  a  constant  function  of  magnitude  1,  while 
x(t)  transforms  to  X(f).  Their  product  is  identically  X(f),  and 
the  inverse  Fourier  -  rar.s  form  yields  x(t).  It  is  obvious  that  any 
input  function  ,  u(t),  will  produce  an  output,  x(t),  identical  to  itself, 
as  it  should,  if  the  ^3:^11  is  transparent.  This  example  should 
not  be  considered  complete  proof  of  the  transfer  function 
approach,  but  merely  demonstrates  an  extreme  situation. 


(if  h(t)  is  the  unit  impulse 


Figure  A.  GRAPHICAL  EXAMPLE  OF  TRANSFER  FUNCTION  APPROACH  AND  "PROOF1 


The  program  developed  to  apply  the  transfer  function  was 
called  XFORM.  The  Cooley-Tukey  fast  Fourier  transform  routine  was 
used  as  the  core  of  this  program.  No  tapers  were  applied  to  the 
time  domain  truncation  function  in  order  to  reduce  leakage  because 
the  frequency  distribution  of  the  waveform  was  unknown.  All  high 
frequency  information  was  desirable,  and  a  taper  could  have 
destroyed  that  information.  Also,  Hanning  or  cosine  windows  often 
smooth  the  waveform.  This  would  art ■ f icially  widen  the  plume;  an 
undesirable  effect.  To  keep  computations  to  a  minimum,  the  number 
of  points  in  the  discrete  Fourier  transform  should  be  small.  Crow 
and  Tewscher  (1983)  determined  the  proper  number  to  be  18,  based 
on  the  instrument  high  frequency  cut-off  and  the  approximate 
airspeed.  The  program  XFORM  therefore  averaged  each  pass  into  an 
18  point  series  before  applying  the  transform.  Each  point, 
therefore,  represented  upwards  of  100  samples.  If  the  measurement 
variability  between  samples  is  considered  independent,  this  would 
decrease  the  statistical  significance  of  measurement  errors 
tenfold.  Considering  the  nature  of  the  noise  (instrument  noise, 
intake  airflow  dynamics,  etc.)  and  the  errors  produced,  the  data 
density  achieved  in  this  experiment  appears  to  be  more  than 
necessary  to  achieve  sufficiently  small  error.  Ten  to  twenty 
samples  per  data  point  would  have  produced  accuracy  to  within  50 
ppt,  an  acceptable  level.  The  18  point  series  was  designed  so 
that  the  records  start  and  end  at  zero  concentration,  with  all 
other  points  non-zero,  to  avoid  introducing  false  high-frequency 
components  due  to  discontinuity  or  background  noise.  The 


untransformed  data  set  was  stored  on  mass  storage  for  comparison 
to  the  transformed  data 

The  program  XFORM  :i«.  ::t  entered  the  transfer  fur.cn  -  ifcrcatir.e. 

The  first  task  in  the  subroutine  was  to  determine  the  transfer 
function.  This  was  accomplished  by  first  transforming  the 
experimental  time  series;  a  simulated  "unit  impulse"  as  the  input 
waveform,  and  the  resultant  measured  SF&  concentration  as  the 
output.  Next,  each  frequency  component's  contribution  to  the 
transfer  function  was  determined  by  dividing  the  input  by  the 
output.  As  implied  in  the  earlier  discussion,  using  an  impulse  as 
input  produces  a  smooth  function  in  frequency  space,  contributing 
information  to  the  transfer  function  from  all  frequency 
components.  Since  the  results  of  the  Fourier  transforms  are 
imaginery  numbers,  their  quotient  is  also  imagery,  as  follows: 

.  Y(f)  (at  +  b  ^i)  a^  +  b1+b2  (b  ^  -  aLb2) 

H(f )  *  -  *  -  *  -  +  -  i  (7) 

X(f)  (a2  +  b2i)  a  ^  a$  +  bg 

is  the  inverse  transfer  function; 

is  the  transform  of  the  laboratory  "impulse" ; 

is  the  transform  of  the  laboratory  output; 

are  the  real  parts  of  the  input  and  output 
transforms ; 


where  H(f)-^ 
Y(f) 
Z(f) 

a  i»a2 


are  the  imaginary  parts  of  the  input  and  output 
transforms . 


The  experimental  time  series  (measured  output)  is  then 
transformed,  and  multiplied  by  the  inverse  transfer  function  to 
yield  the  input  waveform, 


U(f)  »  X(f)H(f) 


-1  a  la2a3+^>l^)2a3  - 

a  5  +  b£ 


bla2b3~a  lb2b3 
a  £  +  b£ 


a  xa  2b  3+b  Lb  2b  3 


a  5  +  b  £ 


b  la  2a  3"a  lb2a3 
a  \  +  b£ 


(3) 


where  U(f)  is  the  tranform  of  "true"  input  waveform? 

X(f)  is  the  transform  of  measured  output  waveform? 
a 3  is  the  real  part  of  the  output  transform? 
b3  is  the  imaginary  part  of  the  output  transform. 

Finally,  an  inverse  transform  yields  the  "true"  input  time 
series . 

XFORM  next  called  the  DELAY  and  ROTATE  subroutines.  These 
subroutines  operated  on  the  coordinates  of  the  pass?  therefore, 
when  no  navigation  information  was  available,  they  were  not  used. 
The  DELAY  subroutine  applied  a  constant  time  delay,  translated  as 
a  shift  in  the  coordinates,  to  account  for  the  lag  time  created  by 
system  dynamics.  This  lag  time  was  obtained  daily  in  situ  tests. 
The  ROTATE  subroutine  corrected  the  concentrations  to  produce 
values  appropriate  to  a  pass  perpendicular  to  the  mean  wind 
direction.  In  almost  all  cases,  this  correction  proved  to  be 
negligible,  since  the  flight  paths  were  usually  within  5  degrees 


of  the  desired  direction. 


The  two  resultant  mass  storage  data  sets  were  called  A1R4 
(untransformed  passes)  and  AIRS  (transformed  passes).  Figures  5 
and  6  show  two  examples  of  untransformed  and  transformed  data. 

The  abscissa  represents  distance  from  first  detection  of  SF6.  The 
apparent  change  in  the  peak  location  upon  transformation  results 
from  the  inherent  time  shift  due  to  the  "smoothing"  of  the  input 
waveform.  In  Figure  5,  the  plume  has  been  significantly  narrowed, 
and  the  mass  conserved  with  an  increase  in  peak  concentration. 
Also,  a  second  mode  appears  which  corresponds  to  a  slight 
inflection  in  the  untransformed  data.  This  demonstrates  the 
usefulness  of  the  transfer  function  approach  for  retrieving  high- 
frequency  information.  Figure  6  displays  a  much  broader  plume 
than  in  Figure  5.  The  transformation  does  not  significantly 
change  this  waveform  shape.  Evidently,  the  transformed  plume 
contains  significant  terms  only  at  frequencies  below  those 
affected  by  the  transfer  function.  Also  note  that  the  peak  shift 
remains,  since  time  shifting  translates  to  phase  shifting  in 
frequency  space,  affecting  all  frequencies. 
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KTOHRR  ^  EXAMPLE  OF  ANALYSIS  STEF  ?  OUTPUT  (NARROW  FLUME ) 


date  :  6/28  time  :  1415:49  pass  27 


FTOtrRE  6.  EXAMPLE  OF  ANALYSTS  STEP  ?  OUTPUT  (WIDE  TT.MME) 


-  Missing  Mini-Ranger  Data 


As  mentioned  in  step  one,  at  times  during  the  experiment  one 
or  both  mini-ranger  signals  were  not  available.  Anticipating  this 
problem,  various  reference  points  were  selected  as  starting  points 
for  plume  transects,  and  the  time  of  intersection  logged  on  the 
SF6  analyzer  strip  chart  by  an  onboard  technician.  With  this 
information,  the  flight  heading,  and  airspeed,  the  plume  coordi¬ 
nates  could  be  estimated. 

The  program  MINIFIX  was  written  to  do  the  necessary  analysis. 
The  heart  of  this  program  is  a  look-up  table  that  lists  passes 
with  missing  coordinates,  and  their  associated  reference  point  to 
plume  peak  distances.  These  distances  were  manually  extracted 
from  the  SF6  analyzer  strip  charts.  Another  table  identifies  the 
reference  points  for  the  various  passes.  These  reference  points 
are  shown  in  Figure  2  and  their  positions  listed  in  Table  3. 

Since  the  reference  point  passage  was  logged  instantaneously 
and  the  strip  chart  data  was  output  by  the  analyzer,  the  inherent 
lag  due  to  the  system  dynamics  (mentioned  in  the  previous  section) 
had  to  be  removed  by  MINIFIX.  Care  was  taken  to  operate  on  the 
untransformed  plume  when  applying  this  correction,  to  avoid  the 
time-shift  due  to  the  transformation  process. 

As  in  the  previous  section,  the  plume  was  rotated  to  a  wind 
direction  perpendicular  alignment.  The  untransformed  and  trans¬ 
formed  data  sets  were  output  to  mass  storage  files  AIR6  and  AIR7 
respectively.  Table  4  lists  all  SF6  cross-sections  archived  at 


NPS,  and  also  identifies  those  passes  with  missing  mini-ranger 
data.  An  example  of  the  data  sets,  AIR4  -  AIR7  (identical 
formats)  is  presented  in  Figure  7. 
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Figure  7. 

EXAMPLE  OF  ANALYSIS  STEPS  2,3  OUTPUTS 
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line  1:  time  in  PST,  pass  numbers  start  at  the  beginning 
of  each  experimental  day 
DWD  is  downwind  distance  from  the  source 
data  -each  profile  has  been  time  center  averaged 
to  18  points 

-SF6  is  concentration  in  PPT 
-all  lengths  are  in  meters 
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Step  4  -  Multi-modal  Gaussian  Fits 


In  plume  dispersion  modelling,  mass  distributions  are  most 
often  described  by  the  familiar  Gaussian,  or  normal,  distribution 
on  the  horizontal  plane.  To  parameterize  these  models,  then,  the 
measured  plume  cross-sections  must  also  be  approximated  in  a 
similar  fashion. 

The  success  with  which  a  Gaussian  shape  approximates  the 
actual  measured  cross-wind  profiles  will,  of  course,  vary  a  great 
deal.  Cross-sections  were  often  skewed,  multi-modal,  or 
" square- shaped " .  This  analysis  step  started  by  determining  the 
standard  deviation  of  the  mass  from  the  mean  position.  In 
discrete  form. 


C«  -  E  f(x.)(x. —x  j ) 2  (10) 

T  i»l 

1  N 

B  *  —  Z  f(x.)(x.-xl)  (11) 

T  i*l 

N 

T  *  E  f(x.)  (12) 

i-1  1 

where  N  is  total  number  of  points; 

f(x^)  is  mass,  or  concentration,  at  the  i  point; 

x.  is  the  cross-wind  position  of  the  i  point; 


These  calculations  were  performed  on  both  the  transformed  and 


untransformed  profiles  of  the  previous  step.  As  expected,  the 
transformed  width  was  always  smaller  than  the  untransformed  value, 
due  to  the  "peak  sharpening"  effect  of  the  transfer  function. 

The  next  task  performed  in  this  analysis  was  a  numerical 
curve  fit  to  the  multi-modal  Gaussian  model,  defined  as  follows: 

n  f“(y-ui)2l 

f(y)  *  E  P.exp  - (13) 

i-1  1  l  2o.2  J 

where  n  is  the  total  number  of  modes; 

ik  is  the  cross-wind  position  of  the  i  mode; 

f(y)  is  the  model  value  at  position  x; 

is  the  peak  concentration  of  the  i  mode; 

o^2  is  the  variance  of  the  i  mode. 

Because  no  unique  solution  to  this  curve  fitting  exists,  the 
program  required  interactive  decisions  for  each  profile.  The  user 
initially  decided  on  the  number  and  cross-wind  positions  of  the 
modes.  The  program  then  selected  the  concentrations  at  those 
positions  to  be  the  model's  amplitude  parameters,  and  calculated 
the  mode  variances  necessary  to  minimize  the  squared  deviations 
from  the  fit.  The  fit  and  observed  profile  were  then  graphically 
displayed  for  the  user.  At  this  point  the  user  could  either 
accept  the  fit,  or  alter  his/her  initial  parameters  to  achieve  a 
more  realistic  model.  Once  satisfied,  the  user  "graded"  the 
profiles  subjectively  in  three  categories:  skewness,  ripple,  and 
overall  goodness  of  fit.  Skewness  refers  to  the  assymetry  of 


the  individual  waveforms  associated  with  each  mode.  Ripple  is  the 
high  frequency  "noise"  introduced  to  the  profiles  through  the  fast 
Fourier  transforms.  Goodness  of  fit  judges  how  well  the  Gaussian 
model  approximates  the  observed  profile.  Table  5  lists  the 
complete  set  of  profiles  and  grades  in  each  category. 

The  results  from  this  analysis  step  are  stored  in  the  mass 
storage  data  set  AIR8.  An  example  of  AIR8  is  supplied  in  Figure  8. 
Examples  of  observed  profiles  and  their  associated  analytical 
forms  is  shown  in  Figures  9-12.  Figure  9  demonstrates  a 
reasonably  well-behaved  profile.  Figure  10  shows  a  bimodal 
distribution.  The  Figures  10  and  11  data  hint  at  an  additional 
mode  in  the  distribution;  however,  the  programmer  decided  to 
ignore  the  minor  peak.  Some  subjectivity  was  inevitable  in  this 
analysis  step.  The  high  frequency  components  in  Figure  12, 
however,  are  ripple,  produced  in  the  FFT.  In  this  case,  the  model 
profile  is  probably  closer  to  reality  than  the  transformed  data. 


Table  5 


Subjectively  Determined  Quality  Analysis 


of  Multi-Modal  Gaussian  Fits 


to  SF6  Cross-Sections 


Grading  System 


GRADE 


skewness  goodness  of  fit  aligned 


Pass  #  ripple  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3 


************************************************************* 
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******** 


skewness 


Pass  #  ripple  peak  1  peak  2 


goodness  of  fit  aligned 


peak  3  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3 


DAY  3  (6-24-82) 


skewness 
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c 

— 

— 

c 

— 

— 

Y 

— 

— 

2 

B 

B 

— 

— 

A 

— 

— 

N 

— 

— 

3 

C 

C 

c 

— 

B 

B 

— 

Y 

Y 

— 

6 

A 

A 

u 

U 

A 

A 

A 

Y 

Y 

Y 

7 

A 

C 

— 

— 

B 

— 

— 

Y 

— 

— 

8 

C 

c 

— 

— 

C 

— 

— 

Y 

— 

— 

9 

A 

B 

— 

— 

B 

— 

— 

N 

— 

— 

10 

A 

B 

— 

— 

C 

— 

— 

N 

— 

— 

12 

A 

B 

— 

— 

B 

— 

— 

N 

— 

— 

13 

A 

a 

u 

— 

C 

A 

— 

Y 

N 

— 

14 

A 

c 

— 

— 

C 

— 

— 

N 

— 

— 

15 

B 

c 

— 

— 

A 

— 

— 

N 

— 

— 

16 

A 

A 

— 

— 

A 

— 

— 

Y 

— 

— 

17 

A 

A 

A 

— 

A 

A 

— 

Y 

Y 

— 

18 

B 

A 

C 

— 

A 

A 

— 

Y 

N 

— 

19 

A 

B 

— 

— 

C 

— 

— 

N 

— 

— 

21 

B 

B 

^  — 

— 

B 

— 

— 

Y 

— 

— 

22 

A 

U 

u 

— 

A 

B 

— 

Y 

N 

— 

23 

A 

B 

— 

C 

— 

— 

Y 

— 

~ 

24************«**+******************* 

************** 

******************************** 

25 

A 

U 

u 

— 

A 

A 

— 

N 

Y 

— 

26 

C 

A 

— 

— 

C 

— 

— 

Y 

— 

— 

28 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

29 

B 

B 

— 

— 

B 

— 

— 

Y 

— 

— 

30 

A 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

31 

B 

C 

— 

— 

C 

— 

— 

N 

— 

-- 

32 

A 

C 

c 

— 

B 

B 

— 

Y 

Y 

— 

33 

C 

A 

— 

— 

C 

— 

— 

Y 

— 

— 

36 

A 

B 

B 

— 

A 

A 

— 

Y 

N 

— 

37 

A 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

38 

A 

C 

c 

— 

B 

A 

— 

Y 

Y 

— 

39 

A 

A 

c 

— 

A 

A 

— 

Y 

Y 

— 

41 

A 

C 

— 

— 

A 

— 

— 

Y 

-- 

— 

42***** 

►  *** 

***********'( 

r****«i 

r ****  *** 

43 

A 

C 

A 

— 

A 

A 

— 

Y 

Y 

— 

45 

A 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

46 

B 

A 

c 

— 

A 

B 

— 

Y 

Y 

— 

47 

A 

B 

B 

— 

B 

A 

— 

Y 

Y 

— 

48 

A 

A 

A 

— 

A 

A 

— 

Y 

Y 

— 

49 

A 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

53 

A 

? 

? 

— 

A 

A 

— 

Y 

Y 

— 

54 

A 

A 

C 

— 

A 

A 

— 

Y 

Y 

— 

55 

A 

C 

— 

— 

A 

— 

— 

Y 

— 

— 

56 

A 

B 

B 

— 

A 

A 

— 

Y 

y 

— 

57 

A 

B 

A 

— 

A 

A 

— 

Y 

Y 

— 

59 

A 

C 

B 

— 

C 

B 

— 

Y 

Y 

— 
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C 

C 

C 

B 

35 

C 

N 

Y 
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DAY  4  (6-25-82) 


skewness 

qoodness  of  fit 

aligned 

Pass  # 

ripple 

peak  ' 

1  peak  2 

peak  3 

peak  1 

peak  2 

peak  3 

peak 

1  peak  2 

peak  3 

1 

A 

B 

A 

Y 

2 

B 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

3 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

4 

C 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

6 

A 

A 

— 

— 

A 

— 

— 

Y 

— 

7 

A 

B 

— 

— 

A 

— 

— 

N 

— 

— 

8 

A 

C 

— 

— 

A 

— 

— 

Y 

— 

— 

9 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

10 

A 

C 

— 

— 

A 

— 

— 

N 

— 

— 

11 

A 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

12 

B 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

13 

A 

U 

u 

— 

A 

A 

— 

Y 

Y 

— 

14 

A 

c 

— 

— 

C 

— 

— 

N 

— 

— 

15 

^g*** 

A 

******** 

c 

******4 

r**  ****★**■< 

******** 

c 

******** 

********4 

r*  ****** 

N 

******4 

r********* 

****** 

17 

A 

c 

c 

c 

c 

Y 

19 

A 

B 

— 

— 

A 

— 

— 

N 

— 

— 

21 

B 

c 

— 

— 

C 

— 

— 

Y 

— 

— 

22 

A 

c 

— 

— 

B 

— 

— 

N 

— 

— 

23 

A 

c 

— 

— 

C 

— 

— 

Y 

— 

— 

24 

A 

B 

— 

— 

A 

— 

— 

N 

— 

— 

25 

C 

C 

— 

— 

C 

— 

— 

N 

— 

— 

26 

C 

C 

— 

— 

A 

— 

— 

N 

— 

— 

27 

B 

C 

c 

— 

A 

A 

— 

N 

N 

— 

28 

A 

U 

u 

u 

A 

A 

A 

Y 

Y 

Y 

29 

A 

U 

u 

— 

A 

A 

— 

Y 

Y 

— 

30 

A 

U 

u 

u 

A 

A 

A 

Y 

Y 

Y 

31 

A 

C 

c 

— 

C 

A 

— 

Y 

N 

— 

32 

A 

c 

A 

N 

34 

B 

u 

u 

u 

A 

A 

A 

Y 

Y 

Y 

35 

A 

— 

— 

C 

— 

— 

Y 

— 

— 

36 

A 

B 

— 

— 

B 

— 

— 

N 

— 

— 

37 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

38 

A 

A 

— 

— 

A 

— 

— 

Y 

— 

— 

39 

A 

U 

u 

— 

A 

A 

— 

Y 

N 

— 

41 

A 

U 

u 

— 

C 

A 

— 

Y 

Y 

— 

42 

A 

U 

u 

— 

A 

A 

— 

Y 

N 

— 

43 

A 

B 

— 

— 

A 

— 

— 

N 

— 

— 

44 

B 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

45 

B 

B 

— 

— 

C 

— 

— 

Y 

— 

46 

A 

C 

— 

— 

A 

— 

— 

Y 

— 

— 

48 

A 

U 

u 

— 

A 

A 

— 

Y 

N 

— 

49 

A 

U 

u 

— 

A 

A 

— 

Y 

Y 

— 

50 

A 

U 

u 

— 

B 

B 

— 

Y 

Y 

— 

51 

A 

A 

— 

— 

A 

— 

— 

Y 

— 

— 

52 

A 

U 

u 

— 

A 

A 

— 

Y 

Y 

— 
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DAY  4  (6-25-82) 

( cont’ d) 


skewness 


goodness  of  fit 


aligned 


Pass  #  ripple  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3 


53  A  C  C  —  C  A  —  Y  Y 

54  A  U  U  —  A  A  —  Y  Y 

55  A  B  —  —  A  —  —  N 

Sg^************************************** ***************************** ************* 
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DAY  5  (6-27-82) 


skewness  goodness  of  fit  aligned 


Pass  #  ripple  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3  peak  1  peak  2  peak  3 


43********************************************************************************* 


44* ********************* ******************* ************************** ************** 
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DAY  6  (6-28-82) 


skewness 

goodness  of  fit 

aligned 

Pass  #  ripple  peak  1  peak  2  peak  3 

peak  t  peak  2  peak  3 

peak  1  peak  2  peak  3 

1 

A 

C 

— 

— 

c 

— 

— 

N 

— 

— 

3 

B 

B 

— 

— 

A 

— 

— 

N 

— 

— 

4 

A 

C 

— 

— 

B 

— 

— 

N 

— 

— 

6 

C 

C 

c 

— 

A 

A 

— 

Y 

Y 

— 

9 

B 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

11 

A 

C 

— 

— 

A 

— 

— 

N 

— 

— 

13 

A 

C 

— 

— 

A 

— 

— 

Y 

— 

— 

15 

A 

C 

c 

c 

A 

A 

A 

Y 

Y 

Y 

17*** 

******* 

18 

A 

c 

— 

— 

A 

— 

— 

Y 

— 

— 

19 

A 

c 

— 

— 

A 

— 

— 

N 

— 

— 

21 

A 

c 

— 

— 

A 

-- 

— 

Y 

— 

— 

22*** 

******* 

******** 

r******i 

2 3* ******************************************************************* ★♦*★*★**♦★**1 

26 

B 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

27 

A 

C 

— 

— 

B 

— 

— 

Y 

— 

— 

29 

A 

A 

A 

— 

A 

A 

— 

Y 

Y 

— 

30*** 

*#♦****4 

********* 

********* 

******* 

******** 

r******i 

33 

C 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

34 

A 

A 

— 

— 

A 

— 

— 

Y 

— 

— 

35 

A 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

36 

A 

C 

— 

— 

A 

— 

— 

N 

— 

— 

37 

B 

C 

C 

— 

B 

B 

— 

Y 

Y 

— 

44 

A 

C 

A 

— 

A 

A 

— 

Y 

N 

— 

46 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

48 

A 

C 

— 

— 

A 

— 

— 

N 

— 

— 

53 

A 

C 

— 

— 

B 

— 

— 

N 

— 

— 

54 

A 

C 

— 

— 

B 

— 

— 

N 

— 

— 

55 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

59 

A 

B 

— 

— 

C 

— 

— 

Y 

— 

— 

60 

A 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

61 

A 

B 

C 

— 

A 

A 

— 

Y 

Y 

— 

62 

B 

C 

— 

— 

C 

— 

— 

Y 

— 

— 

65 

B 

C 

— 

— 

B 

— 

— 

Y 

— 

— 

66 

B 

B 

— 

— 

A 

— 

— 

Y 

— 

— 

67 

A 

C 

— 

— 

A 

— 

— 

Y 

— 

— 

68 

A 

C 

c 

— 

A 

A 

— 

Y 

Y 

— 

71 

A 

— 

B 

B 

— 

N 

N 

— 

72 

A 

C 

— 

— 

C 

— 

— 

N 

— 

— 

73 

A 

A 

A 

A 

N 

Y 

Y 

74 

A 

B 

— 

— 

A 

— 

— 

Y 

- 

— 

77 

C 

C 

c 

c,c+ 

A 

A 

A,A+ 

Y 

Y 

Y,’ 

78 

A 

C 

— 

— 

A 

-- 

— 

Y 

-- 
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date  :  6/25  time  :  1217:11  pass  11 


FIGURE  9  ANALYSTS  STEP  OUTTUT 


date  :  6/22  time  :  1536:22  pass  19 


date  :  6/28  time  :  1415:49  pass  27 


FIGURE  13.  ANALYSIS  STEP  OUTPUT 


date  :  6/25  time  :  1147:30  pass 
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Step  5  -  Calculation  of  Hourly  Averages 


Many  dispersion  models  attempt  to  predict  concentrations 
expected  when  averaged  over  a  one-hour  period.  In  order  to  relate 
the  results  of  this  data  set  to  those  of  the  past,  and  also  in 
order  to  satisfy  contractual  agreements,  this  analysis  step  formed 
one-hour  averages  of  the  horizontal  and  vertical  plume  growth 
parameters,  Oy  and  az.  In  addition,  this  step  added  a  header 
to  the  data  set  containing  a  variety  of  averaged  meteorlogical 
quantities. 

The  basic  procedure  in  this  step  was  to  read  in  half-hour 
average  met  data  twice,  form  one-hour  average  met  data,  read  in 
tracer  data  for  the  current  hour,  bin  the  tracer  profiles 
according  to  range  from  the  release  point,  and  perform  the 
averaging  calculations  for  the  plume. 

The  meteorological  data  was  described  in  Schacher  et  al . 
(1982)  and  was  exclusively  collected  at  the  release  platform.  To 
account  for  plume  flight-time  from  the  platform,  a  lag  of  one-half 
hour  was  applied  when  synchronizing  the  two-data  sets.  Even  with 
this  adjustment,  many  problems  exist  in  determining  the  appro¬ 
priate  meteorology.  Due  to  spatial  inhomogeneity,  meteorological 
conditions  at  the  platform  become  less  representative  of  the 
average  met  conditions  experienced  by  the  plume  as  the  downwind 
distance  to  the  aircraft  transects  increase.  Also,  meteorological 
averages  tend  to  differ  significantly  from  hour  to  hour,  implying 
that  stationarity  through  the  one-hour  period  may  be  a  weak 
assumption. 


For  each  experimentation  day,  four  range  bins  were  selected, 
based  on  the  distribution  of  individual  transect  downwind 
distances.  An  attempt  was  made  to  maximize  the  number  of  passes 
in  each  range  bin  for  all  hours,  while  minimizing  the  standard 
deviation  of  the  downwind  distances  within  each  range  bin.  Table 
6  lists  the  range  bins  for  each  day. 

Table  6 

Range  Bins  for  Hourly  Averages  of  Plume  Parameters 


1 

1 

Transect 

Downwind  Distance 

(m) 

DAY 

I  BIN  1 

BIN  2 

BIN  3 

BIN  4 

6-21-82 

|  0-1000 

1001-2000 

2001-3000 

3001-4500 

6-22-82 

|  0-1000 

1001-2000 

2001-3000 

3001-4000 

6-24-82 

|  0-2000 

2001-3000 

3001-4000 

4001-5000 

6-25-82 

|  0-2000 

2001-3500 

3501-5000 

unused 

6-27-82 

|  0-2000 

2001-3500 

3501-5000 

5001-6500 

6-28-82 

|  0-2500 

2501-5000 

5001-7500 

7501-10,000 

6-29-82 

)  0-2500 

2501-5000 

5001-7500 

7501-10,000 

One 

major  problem 

in  this  analysis  step  was 

collecting  a  suf 

ficiently  large  number  of  transects  for  a  range  bin  during  a  given 
hour.  Typically,  this  number  was  5  to  12  passes  per  average. 


Discussion  of  the  possibilities  and  consequences  of  insufficient  sampling 
will  be  presented  in  a  later  section. 

The  first  averaging  calculations  performed  for  each  hour-range  bin 
were  the  average  and  standard  deviation  of  the  bin's  downwind  distance. 

The  downwind  distance  (DUD)  of  a  cross-section  was  interpreted  as  the  straight 
line  distance  from  the  release  platform  to  the  plume  center.  As  stated  above, 
the  standard  deviations  of  the  DWDs  for  a  range  bin  was  minimized  to  determined 
range  bin  boundaries.  All  DUD  standard  deviations  are  less  than  200  m. 

Five  different  horizontal  plume  parameters  were  calculated  for  each 
hour-range  bin.  Each  operated  on  the  ensemble  of  transects  for  a  bin  in  a 
different  way.  Table  7  gives  symbolic  definitions  used  in  the  discussion 

that  follows. 


Table  7 


Definitions  for  Horizontal  Plume  Parameters 
"  Calculated  tor  Each  Hour/Ranqe  Bin 


Symbol 

Definition 

°yd 

Mean  total  standard  deviations  of 
the  horizontal  mass  distributions 
from  direct  calculations. 

°yf 

Mean  total  standard  deviations  from 
the  uni-modal  Gaussian  fits. 

ayw 

Mean  total  standard  deviations  from 
the  uni-modal  Gaussian  fits  weighted 
by  the  peak  concentration. 

°yt 

Mean  total  standard  standard 
deviations  from  the  uni-modal 
Gaussian  fits  averaged  in  a  fixed 
cross-wind  coordinate  system. 

aym 

Mean  standard  deviations  from  the 
multi-modal  Gaussians  fits. 
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was  the  mean  of  the  standard  deviations  of  the  horizontal 
mass  distributions  as  defined  in  Equations  9-12,  operated  on  the 
transformed  data.  The  cross-wind  coordinate  system  was  allowed  to 
float  in  this  average.  In  other  words,  this  average  is  not 
affected  by  plume  drift. 

Oyf  is  the  analytical  equivalent  of  the  above.  The 
parameters  obtained  during  the  multi-modal  Gaussian  fits 
previous  section  were  combined  to  form  a  single  mode  fit, 
values  averaged.  The  derivation  follows.  In  continuous 
mean  position  of  the  mass  can  be  defined  as  the  expected 
Y,  the  "random  variable"  composed  of  all  y  values. 

e(y)  *  yf(y)dy 

where  Y  is  the  "random  variable"; 

E(Y)  is  the  expected  value  of  Y; 
y  is  the  cross-wind  position; 

F(y)  is  the  density  function  of  y. 

The  variance  is  simply  the  second  moment  of  the  distribution, 
taken  about  the  mean. 

a2  =  E[(Y-u  ?]  =  E(Y2)  -  y2  (15) 

where  u  is  the  distribution  mean;  identically  E(Y) . 

In  the  case  of  the  multi-modal  Gaussian  model,  the  mass,  or 
concentration  distribution,  is  described  by  Equation  13, 
repeated: 


of  the 
and  thos  e 
form,  the 
value  of 

(14) 
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(16) 


n 


f(y)  = 


E  P  .  exp 
i=l 


'(y-lii)  2 


2oi 


where  f(y)  is  the  concentration  at  cross-wind  position  y, 

P^  is  the  peak  concentration  of  the  i  mode. 

The  density  function  can  be  formed  by  simply  normalizing 
Equation  16  by  the  integrated  mass.  The  mean,  or  expected  value, 
of  Y  is  then  easily  derived  as  follows: 


F(y) 


=  f(y)- 


„l*F(y)dy 


n  P  .exp[-(y-u . ) 2/2a . 2] 
E  — = -  1  z 


i=l 


n 

E 

j  =  l 


j2n  a  .  P . 

li 


(17a) 


E(y)  * 


4- os 

/  yF ( y ) dy 


n  a.  P. 
j-1  E  1  1 
i=l 


+•  y  expC-(y-«.)  2/2>  2] 

/  - — — 3 - 3 —  dy  (17b) 

.  J2?  a. 


E(y)  = 


.  i~l  aiPiui 


n 
E 

i=  1 


°iPi 


(17c) 


Rearranging  Equation  15  for  the  ith  mode  yields: 


El(V  )  - 


(18a) 


In  a  similar  fashion  to  the  above  and  using  the  principle  of 
superposition,  it  can  be  shown: 
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E(Y2) 


i-i  aipilui  *  C_V 

"  °ipi 


(18b) 


Again  using  Equation  15,  the  standard  deviation  of  the  ensemble 
profile  with  n  modes  is: 


'  ?  ,  [ o.  P.  (  u?  +  a?  )  -  o .  P.  u  ■  1 

X*1  1  11  1  1  i  i  i 


°iPi 


(19) 


°5£was  obtained  bY  using  Equation  19  for  each  profile,  and 
averaging  all  values  in  each  hour-range  bin.  Results  were  tested 
by  numerically  integrating  the  same  profiles  and  calculating  oy 
as  in  Equations  9-12.  Results  were  within  l*. 

is  weighted  by  the  peak  concentration  of  the  member 

profiles.  This  parameter  is  an  attempt  to  bias  the  mean  value 
toward  the  cloud  width  near  the  plume  centroid  on  the  vertical 
axis,  which  is  ideally  at  the  surface  for  a  surface  release.  If 
o  is  truly  independent  of  height,  o  should  be  indent ical  to 

J  yw 


Oyt is  defined  as  the  mean  total  standard  deviations  from 
the  uni-modal  Gaussian  fits  averaged  in  a  fixed  coordinate  system, 
o  t  was  obtained  in  identical  fashion  to  except  each  transect 
was  fixed  in  the  cross-wind  coordinate  system  before  averaging  so 
that  the  effects  of  plume  centerline  drift  are  included. 


was  consistently  larger  than  u 


The  difference 


o  ^  was  consistent Ay  larger  man  °yf 
between  the  values  can  be  interpreted  as  the  degree  to  which  plume 
meander  dominates  the  hourly  averages.  In  other  words,  a  time- 
averaged  concentration  profile  can  be  divided  into  two  components. 


Plume  spread  due  to  relative  diffusion,  in  which  there  is  no  fixed 
axis,  is  represented  by  mean  fit.  Henceforth,  this  will  be 
called  the  diffusive  component,  and  is  often  referred  to  as  puff 
diffusion.  It  is  chiefly  influenced  by  turbulence  of  length 
scales  close  to  the  size  of  the  cloud.  Plume  spread  due  to  single 
particle  diffusion  relative  to  a  fixed  axis  is  theoretically 
approached  in  Taylor's  (1921)  theorem.  Plume  growth  under  this 
theory  is  influenced  by  the  integrated  energy  spectrum,  or 
turbulence  of  all  scales.  a  t  is  representative  of  this  time- 
averaged  quantity.  The  difference  between  and  oyf  fit  is  the 
time-averaged  plume  spread  due  to  turbulence  of  scales  either  much 
larger,  or  much  smaller,  than  the  cloud  size.  The  later  contribu¬ 
tions  are  negligible.  The  former  turbulence  scales  tend  to  move 
the  whole  instantaneous  plume  in  a  "snake-like"  fashion  and  will 
hereafter  be  referred  to  as  the  meander  component. 

The  final  horizontal  plume  parameter  calculated  was  0  ,  the 

mean  standard  deviations  from  the  multi-modal  fits.  This  quantity 
was  the  mean  of  all  the  individual  mode  widths  in  a  floating  coor¬ 
dinate  system.  The  origin  of  multiple  modes  in  an  instantaneous 
profile  is  yet  unexplained;  therefore,  the  significance  of  this 
calculation  is  unknown.  This  parameter  increases  only  slowly  with 
range,  and  may,  in  fact,  define  the  size  of  coherent  turbulent 


structures 


A  parameter  calculated  during  this  analysis  step  closely 
related  to  horizontal  diffusion  was  the  off-axis  position  of  the 
mean  mass.  This  is  the  difference  between  the  actual  position  of 
the  mean  mass  and  that  position  calculated  from  the  mean  wind 
vector.  The  quantity  shows  any  inhomogeniety  in  the  mean  wind 
field,  such  as  a  sea  breeze's  veer  with  decreasing  distance  to  the 
shoreline.  It  also  reveals  meander  produced  by  motions  of  time 
scales  longer  than  the  one-hour  averaging  period. 

The  vertical  standard  deviation  of  the  concentration  is  not 
measured  instantaneously,  and  therefore  must  be  interpreted  from 
the  horizontal  cross-sections  for  each  hour-range  bin.  This  was 
accomplished,  when  possible,  by  calculating  the  cross-wind 
integrated  concentration  of  each  profile,  and  then  performing  a 
single-sided  Gaussian  fit  in  the  vertical  through  the  data 
points . 

The  cross-wind  integrated  concentration  is  calculated  from 
the  fitted  profiles  and  defined  as  follows: 

n 

+•1 

CWICZ  -  _./  i-1  Piexp[-(y-ui)2/2aJi]dy  (20) 

n 

CWICZ  -  i-1  °yiPi 


(21) 


where  CWIC^is  cros3-wind  integrated  concentration 
in  ppt-m  at  a  height  z, 

Oy^  is  the  standard  deviation  of  the  i  mode  in 

a  given  profile, 

P^  is  the  peak  concentration  of  the  i  mode, 

is  the  mean  position  of  the  i  ^  mode. 


The  model  from  which  oz  was  estimated  is : 


22  ] 


CWIC  «  CW1C o  exp  [- 

2  o2 
z 


(22) 


where  oz  is  the  vertical  standard  deviation  of  mass. 

By  linear  regression  of  ln(CWICz)  versus  z2,  oz  becomes  a 
function  of  only  the  slope,  while  CW1CQ  is  a  function  of  the 
intercept  as  follows. 


'z 


(2a) 


(23) 


CWICQ  »  exp(b)  (24) 

where  a  is  the  slope  of  the  ln(CWICz)  vs  z2  line; 
b  is  the  intercept  of  the  line. 

Errors  in  the  proposed  model  presented  in  Equation  22  can  be 
introduced  by  either  a  differing  vertical  shape  of  the  concentra¬ 
tion  profile  or  a  non-negligible  deposition  of  SF6  onto  the  sea 
surface.  The  profile  shape  was  examined  by  visual  inspection  of 
the  ln(CWICz)  vs  z2  plots.  The  scatter  of  the  points  about  the 
regression  line  appeared  to  be  unbiased  in  the  vertical  for  the 
cases  examined,  indicating  that  the  exp(-z2)  model  was  reasonable 


The  possibility  of  mass  loss  was  examined  by  comparing  the 
ground-level  cross-wind  integrated  concentration  predicted  by 
Equation  24  to  the  value  forced  by  the  source  emission  rate. 

The  Gaussian  plume  model  requires: 


where  CWICG*  is  ground-level  cross-wind  integrated 

concentration  predicted  by  the  Gaussian 
plume  model 

Q  is  the  emission  rate,  25  lb.  SF6/hr, 

a  is  the  range-dependent  vertical  plume  parameter, 
z 

u  is  the  mean  wind  speed. 

Figure  13  shows  the  ratio  of  the  two  values  of  ground- level 
cross-wind  integrated  concentration  as  a  function  of  range. 
Ideally,  this  ratio  should  be  1  for  mass  balance.  Most  points  are 
within  a  factor  of  2.  The  points  are  nicely  scattered  about  the 
identity  ratio,  and  there  appears  to  be  no  range  dependence  from 
0-9  1cm. 

Based  on  these  results,  this  analysis  suggests  that  the 
hourly  averaged  oz  values  determined  by  Equation  22  are 

reasonable. 


Step  6  -  Plume  Parameters  as  a  Function 


of  Range  and  Stability  Class 

This  analysis  step  uses  the  hourly  averaged  tracer  and 
meteorological  data  produced  in  Step  5  to  parameterize  range- 
dependent  plume  growth  as  a  function  of  commonly  obtained 
shipboard  meteorological  measurements .  This  step  uses  only  the 
fixed  fit  <jy  in  the  horizontal  plume  growth  parameterization. 

Future  analysis  will  concentrate  on  some  of  the  other  forms  of  the 
horizontal  plume  dimension,  in  order  to  reduce  scatter  and  examine 
the  effects  of  averaging  time. 

This  analysis  attempts  to  classify  the  plume  properties  on  a 
modification  of  the  well-known  Pasquil 1-Gif ford  table.  (See  Gifford 
[1976]).  The  original  scheme  first  estimates  insolation,  based  on 
cloud  cover  and  time  of  day.  Insolation  range  bin  and  mean 
windspeed  then  determine  the  appropriate  stability  class.  The 
scheme  essentially  makes  use  of  the  strong  relationship  between 
insolation  and  buoyancy  production  of  turbulence  over  land,  while 
relating  mean  windspeed  to  mechanical  turbulence. 

This  scheme  is  not  applicable  over  water  because,  first  of 
all,  buoyancy  is  only  weakly  dependent  on  insolation  over  the 
oceans,  due  to  the  large  specific  heat  of  water.  Air-surface 
temperature  differences,  the  primary  factor  in  buoyant  production 
near  the  surface,  are  more  often  the  result  of  advection  of  either 
water  or  air  masses  than  insolation.  Second,  while  mechanical 


mixing  is  still  primarily  a  function  of  mean  windspeed  over  the 
ocean,  the  analytical  form  of  that  relationship  is  quite 
different. 

In  order  to  find  a  common  link  between  dispersion  over  water 
and  land,  the  fundamental  physical  mechanisms  must  be  examined. 

At  a  given  height,  dispersion  is  primarily  a  function  of  zQ,  the 
characteristic  surface  roughness  length;  and  L,  the  Monin-Obukhov 
length,  defined  as  follows: 


where  u*  is  the  friction  velocity, 

Cp  is  the  specific  heat  of  air  at  constant  pressure 
p  is  the  air  density, 

T  is  the  absolute  air  temperature, 
k  is  the  von  {Carman' s  constant, 
g  is  the  acceleration  of  gravity, 

H  is  the  vertica  heat  flux. 

In  a  now-classic  paper  by  Golder  (1972),  these  quantities 
have  been  related  to  the  Pasquill-Gif ford  stability  classes. 
During  the  BLM  experiments,  Schacher  et  al .  (1982b)  measured  the 
variables  necessary  to  compute  zQ  and  L.  Schacher  et  al. (1982a) 
developed  a  modified  Pasquill-Gif ford  classification  (referred  to 
as  NPS  scheme)  by  relating  zQ  and  L  to  routine  meteorological 
measurements,  and  examined  the  behavior  of  aQ,  the  standard 
deviation  of  the  wind  direction,  as  a  function  of  the  NPS  scheme. 
The  analysis  reported  here  extends  this  concept  one  step  further. 


using  the  NPS  scheme  of  determining  stability  class  together  with 
actual  trace  gas  measurements  to  build  a  family  of  curves. 

The  Schacher  scheme  requires  four  routine  meteorological 
measurements  to  define  stability  class:  mean  windspeed,  relative 
humidity,  air  temperature,  and  sea  surface  temperature.  Three 
sets  of  curves,  for  50%,  80%,  and  95%  relative  humidity,  are  used 
to  determine  the  class.  Figure  14  shows  the  result  for  50% 
humidity.  From  the  air-sea  temperature  difference  and  the  mean 
windspeed,  an  appropriate  Pasquill-Gif ford  stability  class  is 
chosen  by  interpolation  between  curves.  The  complete  set  of 
algorithms  is  presented  in  Table  8.  Two  important  points  are, 
first,  under  this  scheme,  stability  classes  A,  F,  and  G  are  not 
represented  and  second,  the  scheme  breaks  down  at  windspeeds  under 
2  m/s. 

At  windspeeds  under  2  m/s,  unless  conditions  are  highly 
stable,  turbulence,  and  therefore  turbulent  diffusion,  becomes 
highly  inhomogeneous  on  a  horizontal  plane.  Defining  a  stability 
class  in  order  to  define  plume  spread  for  a  Gaussian  dispersion 
model  implies  homogeneous,  steady-state  conditions.  Defining 
stability  class  A  over  the  ocean  is  probably  unnecessary,  and  may 
be  inappropriate  because  it  is  unlikely  the  sea  surface  can  supply 
upward  heat  flux  capable  of  supporting  extreme  super-adiabatic 
conditions.  Defining  classes  F  and  G,  on  the  other  hand,  is 
important  for  coastal  regions.  Kristensen  et  al.  (1981)  gives 
many  over-water  examples  where  these  conditions  prevail  for 
extended  periods  of  time.  Discussion  of  this  problem  is  given  in 
Appendix  B. 


utm/sec) 


Table  8 


P«G  Stability  Class  Scheme 
Adapted  to  a  Marine  Boundary  Layer 

U  =  a0  +  aL4T  +  a2AT2  +  a3AT3  +  a1+ATl+ 

where  a0,  a1#  a2,  a4  are  constants; 

U  is  windspeed; 

AT  is  (air  temperature*  -  sea  surface  temperature)  in  'C 
~Relative  Boundary 


Humidity 

Line 

ao 

al 

a2 

a3 

a  4 

50% 

BC 

1.59318 

-0.95150 

-0.09711 

-0.00610 

-0.00014 

CD 

2.36805 

-1.61613 

-0.18965 

-0.01315 

-0.00031 

DE 

-0.55452 

2.65966 

-0.34382 

0.02783 

-0.00087 

80% 

BC 

1.12799 

-1.08521 

-0.11388 

-0.00707 

-0.00016 

CD 

1.21695 

-2.06787 

-0.25450 

-0.01708 

-0.00040 

DE 

0.56149 

2.53558 

-0.35185 

0.03053 

-0.00100 

95% 

BC 

1.18368 

-0.85413 

-0.05274 

-0.00248 

-0.00005 

CD 

1.12545 

-1.79684 

-0.16237 

-0.00869 

-0.00017 

DE 

0.90463 

2.74354 

-0.47268 

-0.04718 

-0.00165 

*  Optimum: 

10  meter 

measurement 

but  any 

surface  . 

layer  value 

is 

acceptable . 


The  basic  model  used  to  parameterize  plume  growth  for  each 


stability  category  was 


a  j  8 


y,z  ref 


y>z  ref 


where  ay  2(x)  is  the  horizontal  or  vertical  standard 
" s “  deviation  of  the  normally  distributed 

mass  at  range  X  ; 


(27a) 


Oy  _  re£  is  a  constant  for  a  given  stability  category 
representing  an  appropriate  a  at  a  range 

Xytz  ref' 


aj3  are  constants  for  a  given  stability  category 

representing  horizontal  or  vertical  plume 
growth. 


For  comparison  with  accepted  overland  models  of  similar  form, 

Xy  z  ref  was  chosen  to  be  100  m.  To  simplify  notation,  Equation 
27a  can  be  expressed  as  follows: 


o(x)  =  bx 


(27b) 


where  a(x)  is  either  cy  or  crz  , 


is  either  oy  ref/(100f  or  CTZ  ref^100^3'* 


b 

c  is  either  a  or  3 . 

The  regression  analysis  was  performed  in  several  different 
fashions  (to  be  discussed)  for  intercomparison,  but  all  were 
designed  to  minimize  the  mean  fractional  error,  defined  as 
follows: 


2 ( P-0) 

MFE  *  -  (28) 

P+0 
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where  P  is  predicted  plume  parameter, 

0  is  observed  (measured)  plume  parameter. 

Using  this  error  analysis,  instead  of  the  usual  mean  square 
error,  gives  logarithmically  unbiased  results;  an  over-estimate  of 
n  x  measured  value  is  the  same  as  an  underestimate  of  1/n  x 
measured  value.  This  implies  that  overpredictions  are  more 
heavily  weighted  than  underpredictions.  This  is  a  desirable 
trait,  since  the  data  set  has  a  lower,  but  no  upper  boundary. 

Also,  the  standard  deviation  of  the  MFE  is  a  measure  of  the  preci¬ 
sion  (scatter)  of  the  estimate;  another  useful  characteristic. 
Irwin  (1982)  gives  a  similar  example  of  the  use  of  MFE  in  a 
sensitivity  analysis  of  overland  dispersion  models. 

Equation  27b  contains  two  unknowns.  The  coefficient  b 
essentially  represents  the  initial  conditions,  or  short-range 
diffusion,  which  has  not  been  measured  directly  over  the  ocean. 

The  exponent,  c,  represents  the  curvature  of  the  scatter  plot,  or 
the  deviation  from  linearity  of  plume  growth  as  a  function  of 
range.  Regressing  ln(a(x))  versus  ln(X)  and  allowing  both  b  and  c 
to  vary  will  not  yield  a  unique  solution.  However,  selecting  a 
discrete  set  of  values  for  either  b  or  c  will  produce  a  single  MFE 
minimum. 

The  first  regression  scheme  attempted  was  to  select  a 
discrete  set  of  values  for  c  and  examine  the  standard  deviation  of 
the  MFE.  In  all  cases  varied  only  slightly,  suggesting  that 

there  was  no  preferred  combination  of  b  and  c. 


Next,  c  was  held  constant  and  b  allowed  to  vary.  The  value 
of  c  was  chosen  to  be  0.85  for  horizontal  diffusion.  This  was 
based  on  a  sensitivity  study  of  various  models  as  they  apply  to 
the  Brookhaven  over-water  oil  smoke  experiment  conducted  off  Long 
Island  (Michael  et  al.  [1973]).  The  study  suggested  the  0.85 
value  to  be  appropriate  for  all  stability  classes.  Over  land,  c 
varies  from  about  0.80  for  Pasquill-Gif ford  (P-G)  class  E  to  1.00 
for  classes  G-A.  In  most  cases,  holding  c  constant  produced 
unreasonable  values  of  b.  In  other  words,  the  model 
misrepresented  short-range  diffusion. 

To  remedy  the  problem,  the  approach  was  reversed;  estimates 
of  short-range  diffusion  were  assumed  and  the  curvature  term 
forced.  As  previously  mentioned,  no  short-range  diffusion  data 
are  available.  However,  statistical  theory  introduced  by  Taylor 
(1921)  and  applied  by  Pasquill  (1971)  and  Draxler  (1976)  allow 
estimates  of  short-range  diffusion.  Specifically,  in  the 
horizontal  case. 


oy(T)  - 


(29) 


where 


is  the  standard  deviation  of  the  cross-wind 
velocity  component; 

is  the  diffusion  time; 

is  a  universal  function; 

is  the  Lagrangian  time  scale. 
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Approximating  ayT  **  a0x 


ay(x)  =  o0xfy  (30 

where  a.  is  the  standard  deviation  of  the  wind  direction, 
x  is  the  downwind  distance. 

Sheih  (1981)  has  experimentally  determined  the  "universal" 
function  over  Lake  Michigan  for  various  P-G  categories  from 
trajectories  of  neutral ly-buoyant  balloons  in  the  surface  layer. 
Sheih  (1981)  used  the  model: 


where  t'  is  an  "apparent"  integral  time  scale. 

Li 

Table  9  lists  Sheih' s  experimentally-determined  "apparent" 
integral  time  scales  and  Draxler's  overland  equivalent.  Draxler 
only  separated  data  into  stable  or  unstable;  therefore,  no  "D" 
value  is  presented.  Notice  the  large  time  scale  in  neutral 
conditions,  representing  a  large  "memory"  of  an  air  parcel's 
trajectory.  This  is  probably  a  response  to  synoptic  scale 
disturbances.  In  non-neutral  conditions,  the  time  scale  is 
significantly  less  than  the  over-land  counterpart. 

Equations  30  and  31  can  be  used  with  measured  values  of  a g 
to  obtain  horizontal  short-range  parameters.  The  a 9  values 


Table  9 


Sheih*s  Apparent  Integral  Time  Scales 


P-G  CLASS 

C _  D  _ E 

(all  values  in  seconds) 


horizontal 

372  £  29(617) 

4056  -  223 

70.3 

3.2(617) 

vertical 

10.6  i  1.1(309) 

31.5  £  2.1 

21.7 

+ 

1.3(617) 

(  )  Draxler's  over-land  results 


Table  10. 

* 

Horizontal  and  Vertical  Wind  Variance  Values  from 


Class 

Central 

California  Air  Quality 

Studies  III 

and  IV 

♦  Hrs 

1  Hr^ 

1  Hr(^ 

1  min  0ft 

Oq  Gifford 

B 

10 

31.0 

11.8 

7.2 

20 

C 

10 

17.3 

9.8 

7.3 

15 

D 

129 

9.1 

3.3 

2.6 

10 

E 

36 

12.6 

1.5 

2.1 

5 

*  all  values  in  degrees,  and  measured  at  10  m. 


(35) 


(36) 


(36) 


fM  -  0 .4  [6 . 78  +  2 . 39  ( z/L  -  2.0)  ]  for  z/L^2.0 


-.333 

o  “  h  for  h  >  333  m 


-.175 

a  =  JLil  h  for  z/L<0,  25<h<333  m 

♦  u 


where  L  is  the  Monin-Obukhov  length; 
u*  is  the  friction  velocity; 
u  is  10  m  windspeed; 

h  is  the  inversion  height. 


As  mentioned  above,  the  reference  distance  used  for  the  short 
range  diffusion  parameter,  z  ref ,  was  100  m.  At  this  range, 
Equations  30-31  produce  the  results  presented  in  Table  11.  The 
minimum  and  maximum  values  result  from  deviations  in  the 
"universal  function"  due  to  uncertainties  in  the  diffusion  time 
(windspeed)  and  the  apparent  integral  time  scale  (error  margins  in 
Table  9).  Sheih  (1981)  did  not  present  a  value  of  t£for  class  B; 
therefore,  values  of  Table  11  are  based  on  "reasonable"  t^ 
values. 


Table  11 


Calculated  o  -  values  at  100  m. 

_  y,z  ref  _ 


oy  ref 


ffz  r 


Class 

min 

mean 

max 

min 

mean 

max 

B* 

21.65 

27.01 

32.48 

6.17 

8.23 

10.29 

C 

24.39 

25.90 

27.10 

6.99 

8.70 

10.23 

D 

14.77 

15.09 

15.41 

3.20 

3.73 

4.19 

E 

14.35 

16.11 

17.44 

1.34 

1.61 

1.83 

*  only 

approximate 

An  interesting  aspect 

of  these 

results 

is  that,  for 

the 

horizontal  case 

i,  the  class 

0  and  E 

cases  are 

very  similar 

.  This 

is  the  result  of  compensating  influences  of  a0  and  t^;  the  smaller 
a  values  in  class  D  are  offset  by  the  larger  integral  time  scale 
(memory) . 

With  the  coefficient  term  of  Equation  28  defined,  the 
exponent  can  be  forced  in  the  regression  analysis  scheme.  The 
results,  the  applications,  and  limitations  are  presented  in  the 
next  chapter. 


i 


CHAPTER  II  -  PRELIMINARY  RESULTS 
Additional  Data  Sets 


Three  additional  data  sets  have  been  convolved  with  the  data 
set  described  in  this  report  (see  Table  12).  All  experiments  were 
conducted  with  continuous  surface  releases  of  the  inert  gas  SF6. 
This  implies  that  the  parameterizations  derived  will  be  most  ap¬ 
plicable  to  a  similar  release.  In  addition,  Dabberdt  et  al.(1983) 
produced  some  shoreline  ay  and  az  values  from  the  fourth 
Central  California  experiment  (BLM  IV)  which  are  also  incorporated 
into  our  data  set.  The  first  Gulf  of  Mexico  experiment  (GULF  I) 
was  conducted  during  the  summer.  The  warm  Gulf  water  produced  the 
only  P-G  class  B  and  C  conditions  that  coincided  with  tracer 
releases.  The  third  Central  California  experiment  (BLM  III)  and 
GULF  II  were  conducted  in  winter.  Cool  evening  temperatures 
produced  some  unstable  conditions  during  BLM  III,  but  these  events 
rarely  coincided  with  tracer  releases.  GULF  II  was  conducted 
during  a  stable,  foggy  period. 


Table  12 


ADDITIONAL  DATA  BASES  FOR  OVERWATER,  MEDIUM- RANGE , 
SURFACE-RELEASE  PLUME  DISPERSION  PARAMETERIZATION 


EXPERIMENT 

DATE 

LOCATION 

REFERENCE 

Gulf  of  Mexico 

Air  Quality  Study  I 

Jul 

81 

Cameron,  LA 
(area) 

Dabberdt  et 
al.  (1982) 

Central  California 

Air  Quality  Study  III 

Dec 

81 

Pismo  Beach,  CA 
( area) 

Dabberdt,  et 
al.  (1983) 

Gulf  of  Mexico  II 

Feb 

82 

Cameron,  LA 
(area) 

Dabberdt,  et 
al  (1982) 

71 


The  complete  set  of  additional  data  and  method  of  measurement 
is  supplied  in  Table  13.  Meteorological  data  is  not  tabulated,  but 
stability  categories  were  obtained  in  the  manner  described  in  this 
report. 

Table  13 

ADDITIONAL  1  HR  AVERAGE  PLUME  PARAMETERS 

-all  values  in  meters 

-all  az  values  from  aircraft  transects 
-"s"  indicates  shoreline  collectors  for  a y 
-"a"  indicates  aircraft  transects  for  <jy 
-"b”  indicate  grab  bag  samplers  from  boat 

Experiment  Method  Date  HR  <*y  35  Range 


Exp.  Method 

Date 

HR 

05 

O* 

Ranqe 

BLM  III  s 

12-17-81 

12 

216.5 

6380 

s 

13 

231 

17.5 

6510 

s 

14 

332 

- 

6380 

8 

15 

677 

116 

6630 

8 

16 

299 

39 

6860 

8 

17 

154 

22.5 

6960 

8 

18 

387 

- 

7390 

8 

12-14-81 

12 

194 

18 

6510 

8 

13 

200 

22.5 

6590 

S 

14 

18  7 

23.5 

6530 

S 

15 

176 

12 

6600 

8 

16 

224 

- 

6740 

8 

17 

784 

- 

7310 

8 

12-15-81 

12 

601 

79.5 

7030 

8 

13 

346 

42 

6930 

8 

15 

723 

14.5 

6560 

8 

16 

268 

16.5 

7010 

8 

17 

458 

35.5 

7430 

S 

18 

812 

8290 
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EXP 

Method 

Date 

HR 

Sh 

% 

Range 

BLM  IV 

s 

6-25-82 

11 

117 

6220 

s 

12 

219 

30 

6220 

s 

13 

260 

55 

6220 

s 

14 

239 

36.5 

6220 

s 

15 

149 

46.5 

6240 

s 

16 

156 

6260 

s 

17 

525 

6430 

s 

6-27-82 

11 

139 

6820 

s 

12 

83 

6610 

s 

13 

131 

6670 

s 

14 

202 

34 

6630 

s 

15 

156 

39 

6650 

s 

16 

17  2 

32 

6720 

s 

17 

263 

32 

6640 

GULF  I 

s 

7-20-81 

13 

55 

58.5 

7019 

b 

13 

483 

8661 

s 

14 

671 

9275 

b 

14 

85 

7480 

s 

15 

2088 

8330 

b 

15 

305 

6209 

EXP 

GULF  II 


Method 

Date 

HR 

5* 

Range 

s 

2-23-82 

10 

498 

- 

7847 

s 

11 

238 

76.7 

7724 

a 

11 

146 

4265 

a 

11 

349 

- 

4360 

a 

11 

139 

- 

4662 

s 

12 

471 

33 

8035 

a 

12 

109 

- 

4631 

s 

12 

145 

4553 

a 

12 

115 

463  3 

s 

13 

179 

53 

7741 

a 

13 

89 

■ - - 

4343 

a 

13 

198 

- - 

4370 

a 

13 

163 

_ 

4411 

s 

14 

117 

7912 

s 

15 

315 

57.7 

7984 

a 

15 

295 

4545 

a 

15 

268 

- 

4546 

a 

15 

490 

— 

4044 

s 

16 

489 

46.7 

7309 

a 

16 

395 

— 

4106 

a 

16 

313 

- 

4105 

a 

16 

490 

— 

4044 

s 

17 

107 

33.5 

7494 

a 

17 

116 

- 

7186 

a 

17 

99 

_ 

7123 

s 

18 

101 

_ 

7505 

Vertical  Dispersion  Parameters 


The  encouraging  results  of  the  vertical  dispersion 
parameterization  are  the  well-behaved  form  of  oz  and  the 
distinct  difference  between  classes  D  and  E.  The  discouraging 
aspect  is  that  this  data  contains  no  class  B  or  C  values  for 
az.  Figure  15  shows  the  BLM  IV  scatter  plots  and  regression 
curves  for  classes  D  and  E.  Numerical  results  are  presented  in 
Table  14.  Also  shown  is  the  Turner  (1970)  overland  curves  for 
comparison.  The  figure  shows  obvious  differences  between  classes 
and  a  general  slower  overwater  growth  compared  to  its  overland 
counterpart.  The  slower  vertical  growth  is  physically  realistic 
when  we  consider  surface  roughness.  Lower  values  of  z0  overwater 
produce  smaller  vertical  velocity  fluctuations  during  stable  and 
neutral  conditions,  and  therefore  smaller  plume  parameters.  The 
additional  data  sets  were  not  included  in  the  regression  analysis . 
The  Gulf  data,  by  the  author's  admission,  showed  serious  mass 
balance  problems.  Both  data  sets  were  based  on  airplane  transects 
over  the  shoreline,  where  the  internal  boundary  layer  could  have 
altered  results.  Nonetheless,  this  data  is  included  in  Figure  16 
for  review,  and  supports  our  results.  As  stated  above,  tracer 
data  did  not  coincide  with  periods  of  class  B  or  C  stability. 
Meteorological  data,  however,  was  logged  for  20  complete  hours 
during  these  conditions  (10  hours  apiece).  Based  on  the 
calculated  vertical  wind  variance  for  these  classes,  and  the 
well-behaved  vertical  dispersion  in  the  neutral  and  stable 


FIGURE  15.  ONE  HOUR  AVERAGE  VERTICAL  PLUME  PARAMETER  FROM  CCAQ  IV 


dashed  line  is  Pasquill-Gif ford 
solid  line  is  table  (13) 
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FIGURE  16.  ONE  HOUR  AVERAGE  VERTICAL  PLUME  PARAMETER  FROM 
CCAQ  IV  AND  TABLE  (11)  DATA  SETS 

dashed  line  is  Pasquill-Gif ford 


solid  line  is  Cable  (13) 
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Table  14 


ONE-HOUR  AVERAGE  PARAMETERIZATION  FOR  OVERWATER,  SURFACE- RELEASE 


25.0 

19.0 

10.0 

11.0 

0.75a 

1.00 

0.75a 

1.0 

20.0 

12.5 

8.0 

7.5 

0.70a 

1.00 

0.70a 

0.90 

15.1 

8.0 

3.2 

4.5 

0.69 

0.90 

0.65 

0.85 

16.1 

6.0 

1.8 

3.5 

0.65 

0.80 

0.62 

0.80 

insufficient  data  for  verification 
moderate-range  is  0.1-12  km 
**  Overland  values  from  OTIC  (1980) 


categories,  the  shape  of  the  <jz  curve  is  postulated  in  Table  14. 
Verification  will  proceed  as  unstable,  overwater  data  become 
available  to  the  NPS  Environmental  Physics  Group. 

Horizontal  Dispersion  Parameters 
The  hourly  averaged  horizontal  tracer  data  for  P-G  classes  D 
and  E  with  regression  lines  are  shown  in  Figure  17.  These  results 
are  aesthetically  less  pleasing  than  the  vertical  case  because  of 
the  increased  scatter,  but  some  differences  between  cases  are 
noteworthy.  First,  the  increased  short-range  diffusion  due  to 
meander  for  class  E,  predicted  by  the  theory  of  the  previous 
section,  appears  to  be  realistic  when  examining  the  clusters  in 
the  1-2  km  range.  Second,  clusters  at  greater  ranges  suggest  the 
overall  larger  diffusion  under  class  D  conditions.  The  difference 
is  small,  however,  and  the  parameterizations  of  Table  13  reflect 
this  fact.  As  with  the  vertical  data,  P-G  classes  B  and  C  were 
insufficiently  dense.  Ten  data  points  were  available  in  class  C, 
seven  in  class  B,  and  all  data  were  from  GULF  I.  No  regression 
was  attempted  on  these  data,  and  the  values  in  Table  13  were 
hypothesized,  based  on  the  meteorological  (a.)  data.  Verification 
is  needed. 


FIGURE  17.  ONE  HOUR  AVERAGE  HORIZONTAL  PLUME  PARAMETERS 
FROM  CCAQ  IV  AND  TABLE  (II)  DATA  SETS 

dashed  line  is  Pasquill-Gif ford 


solid  line  is  table  (13) 


As  with  most  tracer  data,  the  points  were  widely  scattered 
about  chosen  regression  lines.  This  characteristic  feature  can  be 
partially  attributed  to  the  highly  variable  nature  of  turbulence 
in  the  atmosphere.  Another  factor  that  significantly  increases 
scatter  for  horizontal  data  is  the  large  energy  in  the  low 
frequency  part  of  the  horizontal  velocity  spectra.  While  a  formal 
spectral  analysis  of  the  wind  time  series  was  not  performed, 
variance  did  significantly  increase  with  longer  sampling  windows , 
up  to  one  hour.  The  time  series  also  suggests  that  this  trend 
would  have  continued  with  a  larger  window.  A  variety  of  overland 
experiments  have  observe  I  lar^e  horizontal  wind  variance  during 
stable  conditions  [Hanna  (1981),  Olesen  et  al.  (1983),  Sagendorf 
and  Dickson  (1974)].  Spectral  analysis  by  Hanna  (1981)  indicated 
a  low  frequency  peak  at  approximately  0.5  hour1'1.  Olesen  et  al. 
(1983)  describe  large  contributions  to  the  energy  spectrum  at 
frequencies  as  low  as  0.35  hour"1.  Kristensen  et  al.  (1981) 
described  increased  plume  meander  in  very  stable  conditions 
resulting  from  these  low  frequency  oscillations,  and  finds  an 
inverse  relationship  with  the  mean  windspeed  (see  Appendix  B). 

Based  on  the  above  references,  it  is  not  surprising  to  find  a 
large  meander  component  in  the  class  E  ay  values.  It  is 
somewhat  unexpected  to  find  a  large  meander  component  in  near¬ 
neutral  (class  D)  stability.  These  findings  are  supported  in  part 
by  Sheih's  (1981)  large  Langrangian  time  scales  in  these  condi¬ 
tions,  which  he  has  suggested  is  the  result  of  "large  scale 


motions 


RD-A148  632  MEASURED  PLUME  DISPERSION  PARAMETERS  OVER  HATER  VOLUME 
1(U)  NAVAL  POSTGRADUATE  SCHOOL  MONTEREV  CA 
C  E  SKUPNIEWICZ  ET  AL.  SEP  84  NPS-61-84-012 


2/2 


UNCLASSIFIED 


F/G  13/2 


NL 


APPENDIX  A 


FORMAT  AIR2 


REDUCE  AIR3 


XFORM  AIR4 


XFORM  AIRS 


raw  data  1 


2 


3-end 


calibrated  1 

data,  rec¬ 
tangular  co-  2 

ordinates 

3 

4 

5 

6 

7 

8-11 

12 

13 

14 -end 


untransforraed,  1 
arranged  data  2 

3 


4 

5 

6-24 


transformed, 
averaged  data 


code,  date  time, 
pass#,  data  quality 
index,  average 
altitude 

code ,  elapsed  time , 

#of  points,  #  of 
transponder  1  pts, 

#  of  transponder 
2  pts,  #  bad  pts 
code,  mini-ranger#, 
mini-ranger  distance, 
analyzer  output 
date, time,  pass# 
#points 

plane  heading,  air¬ 
speed 

wind  direction, wind 
speed 

standard  deviation  of 
output  data 
cross-wind  integrated 
concentration 
transect  altitude 
distance  from  release 
mini-ranger  statistics 
Release  coordinates 

elapsed  time, 
running  plume  width, 
e-w  coordinate,  n-s 
coordinate,  concen¬ 
tration 

date, time,  pass# 
altitude,  windspeed, 
wind  direction 
plane  heading,  air¬ 
speed,  distance 
from  release 
release  coordinates 


elapsed  time,  run¬ 
ning  plume  width,  e- 
w  coordinate,  n-s 
coordinate,  concen¬ 
tration 

Same  as  AIR4 


TABLE  1A  -  (cont'd) 


Output 

Analysis 

Program 

Data  Set 

Brief 

Step 

Name 

Name 

Description 

line# 

Ordered  Content 

3 

MZNZFZX 

AIR6 

corrected 
coordinate s- 
untransformed 
data 

same  as  AXR4 

AZR7 

corrected 

coordinates- 

transformed 

data 

same  as  AIR4 

4 

FIT 

AIRS 

multi-modal 

1 

null 

Gaussian  fit 

2 

plane  heading,  time 

3 

altitude ,  distance 
from  release 

4 

e-w  coordinate , n-s 
coordinate 

5 

width  position 
of  mean  mass 

6 

standard  deviation 
about  mean 

7 

total  plume  width 

8 

null 

9 

null 

10 

peak#l  value,  peak# 2 
value,  etc. 

11 

peak#l  position, peak# 
2  position,  etc 

12 

peak#l  standard  dev. , 
peak# 2  st.  dev.,  etc. 

5 

BOTH 

AZR9 

hourly 

See  Appendix  C  for  complete 

averages 

AZR9  output 

6 

BOTH 

AIR12 

AIR9 

1 

date,  hour,  relative 

condensed 

humidity,  wind 

direction,  sigma 
theta 


2  windspeed,  air  tem¬ 
perature  ( 1 Om. ) ,  sea- 
surface  temperature, 
10/L,  inversion 
height 

3  1st  average 
downwind  distance 
(DWD) ,  1st  standard 
deviation  of  DWD,  1st 
#  of  passes ,  1st  mean 
total  sigma  y 

4  1st  mean  waveform 
sigma  y 
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gr-.L 


aaai 


Brief 

Description 


line* 


AIR9 

condensed 


4  cont.  1st  fixed 

mean  total  sigma  y 
from  fits,  1st  off- 
axis  position  of  mean 
mass,  2nd  average 
DWD,  2nd  st.  dev.  of 

5  DWD 

2nd  #  of  passes,  2nd 
mean  total  sigma  y, 

2nd  mean  waveform 
sigma  y,  2nd  fixed 
mean  total  sigma  y 
from  fits,  2nd  off- 
axis  position  of  mean 

6  mass 

3rd  av.  DWD,  3rd 
standard  deviation  of 
DWD,  3rd  #  of  passes, 
3rd  mean  total  sigma 
y,  3rd  mean  waveform 

7  sigma  y. 

3rd  fixed  mean  total 
sigma  y  from  fits, 

3rd  offaxis  position 
of  mean  mass,  4th 
average  DWD,  4th  st. 

8  dev. of  DWD, 4th#  pas. 
4th  mean  total  sigma 
y,  4th  mean  waveform 
sigma  y,  4th  fixed 
mean  total  sigma  y 
from  fits,  offaxis 
position  of  mean  mass , 
1st  mean  total  sigma 

9  y  from  fits 

1st  weighted  mean 
total  sigma  y  from 
fits,  1st  sigma  z, 

1st  crosswind  inte¬ 
grated  concentration 
(CWIC),  2nd  mean  to¬ 
tal  sigma  y  from  fits 
2nd  weighted  mean 
total  sigma  y  from 

10  fits 

2nd  sigma  z,  2nd  CWIC 
3rd  mean  total  sigma 
y  from  fits,  3rd 
weighted  mean  total 
sigma  y  from  fits, 


Output 

•;  Analysis 

Program 

Data  Set 

Brief 

Marne 

Description 

Ordered  Content 

AIR12 


AIR  9 
condensed 


12,13 


3rd  CWIC,  4th  mean 
total  sigma  y  from 
fits,  4th  weighted 
mean  total  sigma  y 
from  fits,  4th  sig¬ 
ma  z,  4th  CWIC 
null 


Mote: 


Identically  formatted  over-water  data  sets  for  Central  California  Air 
Quality  Exp  III  and  the  two  Gulf  of  Mexico  experiments  are  also  on  file. 
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APPENDIX  B 


OVER-WATER  PLUME  DISPERSION  IN  VERY  STABLE  CONDITIONS 
As  stated  in  the  main  text,  very  stable  conditions  are  not 
uncommon  over  the  ocean.  These  conditions  typically  occur  when 
the  marine  boundary  layer  capping  inversions  lowers  to  the  sea 
surface.  Under  such  conditions,  the  only  true  measure  of 
stability  is  the  atmospheric  temperature  lapse  rate  through  the 
inversion.  Dispersion  in  these  conditions  departs  radically  from 
traditional  turbulent  diffusion  ideas.  Kristensen  et  al.  (1981) 
gives  an  elaborate  theoretical  discussion  of  the  physics  of 
dispersion  in  very  stable  conditions,  identifying  the  key 
parameters  as  averaging  time  and  mean  windspeed.  using  over-water 
tracer  data  at  a  20km  range,  Kristensen  found 

1/3  -4/S 

oy  *  3700  T  U  (1A) 

where  T  is  average  time; 

U  is  mean  windspeed. 

This  formula  is  only  valid  at  20  km,  and  therefore  is  of 
little  value  to  us,  but  demonstrates  the  convincingly  changed 
character  of  diffusion  in  very  stable  conditions. 


CENTRAL  CALIFORNIA  AIR  QUALITY  EXPERIMENT 


Measured  Plane  Dispersion  Parameters  Over  Water:  Volume  2) 


available  on  request  only 
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